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Abstract 

Sgr A East is a supernova remnant located within few parsecs 
from the Galactic Centre (GC). There are good reasons to believe 
that this object is the source of the 7-ray excess detected by HESS 
in the direction of the GC meaning that Sgr A East is likely to be 
an efficient Cosmic Ray (CR) accelerator. Some observations suggest 
that strong magnetic fields may be present in that region which may 
allow the acceleration of composite nuclei in Sgr A East beyond the 
EeV. We show that, if this is case, EeV neutrons should be effectively 
produced by the photo-disintegration of Ultra High Energy nuclei onto 
the infrared photon background (with temperature ~ 40 K) in which 
Sgr A East is embedded. Neutrons with such an energy can reach 
the Earth before decaying and may be detectable under the form of 
a CR point-like excess in the direction of the GC. We determine the 
expected energy spectrum and the amplitude of this signal showing 
that it may be measurable by the AUGER observatory. 



1 Introduction 

One of the main puzzles in cosmic ray physics concerns the origin of cosmic 
rays (CRs) with energies between the knee (-Eknee ~ 3 x 10 15 eV) and the 
ankle (-E an kie ~ 6 x 10 18 eV) of the CR spectrum. It is generally accepted 
that CRs with energy up to the knee are accelerated by the 'so-called' Fermi 
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first order mechanism (for a review see e.g. pQ) in the supersonic outflow of 
galactic SuperNova Remnants (SNRs). This wisdom is motivated not only 
by the well known coincidence between the galactic SN energy output and 
that of galactic CRs but also by the slope and the regularity of the CR energy 
spectrum over several energy decades. The observation of radio and 7-ray 
emissions from a number of SNRs provides further evidence of the validity of 
this scenario. It is more uncertain, however, what is the maximal energy E max 
at which particles can be accelerated in this kind of sources. The continuity 
of the CR spectrum up to ankle suggests that the origin of CR above the 
knee is closely related to that below the knee meaning that -E max may be 
larger than the EeV. 

One of the crucial quantities which enters to determine E mgx is the mag- 
netic field strength in the region where the SNR shell is expanding through. It 
is known that, if the mean value of the galactic magnetic field (-B ga i — 5 /xG) 
is adopted, the theoretically estimated E max cannot exceed the Lagage and 
Cesarsky bound E^ x ~ Z x 10 16 eV 0. Stronger magnetic fields may how- 
ever be found in some dense regions of the Galaxy, like the Galactic Centre 
(GC) [3] or in dense molecular clouds 4j. It is also possible that non-linear 
effects, e.g. cosmic-ray streaming instabilities [Hj, give rise to an amplification 
of the magnetic field in the acceleration region allowing to reach CR energies 
exceeding the Lagage and Cesarsky bound. Even in this case, however, the 
largest CR energies are expected to be reached in SNRs placed in regions 
with the strongest magnetic field backgrounds. 

It would be extremely interesting to find an observational evidence show- 
ing that CR acceleration up to Ultra High Energies (UHEs) (E > 10 18 eV) 
takes place in some galactic SNR. Unfortunately, even EeV protons are sig- 
nificantly deflected in the galactic magnetic field so that the angular cor- 
relation between the arrival direction of these particles and the position of 
their sources should be spoiled. A possible way out may, however, be found 
if EeV secondary neutrons are efficiently produced by the scattering of the 
UHE primary nuclei with the matter and the radiation surrounding the SNR. 
Neutrons with such an energy can travel over galactic distances without un- 
dergoing significant decay and deflections. Although high energy neutrons 
and protons give rise to indistinguishable showers in the atmosphere, a galac- 
tic neutron source may be recognised by the extensive air shower experiments 
under the form of a localised excess of CRs with energy > 1 EeV. Clearly, 
a complete identification of the source requires a matching between the posi- 
tion of the CR excess and that of a radio, X-ray and, possibly, 7-ray emitting 
SNR 0. 

The SNRs which are the most promising sources of EeV neutrons are 
those residing in the densest regions of the Galaxy, where strong magnetic 
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fields as well as thick radiation and gas targets can presumably be found. 
Active SNRs in the proximity of the Galactic Center (GC) and/or in the 
nearby of dense molecular clouds are natural candidates to the role of EeV 
neutron factories. 

In this respect, as it was recently suggested in 0, a particularly interest- 
ing system is Sgr A East. Sgr A East is a SNR located within few parsecs 
from the dynamical centre of the Galaxy. Detailed radio and X-ray obser- 
vations jHl E| allowed to conclude that this object is a mixed-morphology 
(radio shell + X-ray emitting core) remnant of a single type II SN explosion 
which took place between 10 4 and 10 5 years ago (the light propagation time 
is subtracted). This SNR is embedded within a dense ionised gas halo and 
in a background of Infrared (IR) radiation due to dust emission at the tem- 
perature T ~ 40 K pm ITT]. Furthermore, the expanding Sgr A East radio 
shell is interacting with a dense molecular cloud. Strong magnetic fields, as 
large as few milliGauss, have been observed in that region [Hj. 

A 7-ray emission has been also recently observed by the HESS Cerenkov 
telescope in the direction of the GC jT^j. Although the limited angular 
resolution reachable with this kind of observations did not allow a firm iden- 
tification of the source of this emission, several arguments points to Sgr A 
East as the most plausible source. The most convincing among these argu- 
ments is based on the energetic of the emission observed by HESS. In Sec. 0] 
we show, in fact, that the 7-ray flux measured by HESS practically coincides 
with that expected from the decay of neutral pions produced in the hadronic 
collisions of UHE protons which are shock accelerated by a SNR with the 
same characteristics of Sgr A East. On the basis of this argument, as well 
as other arguments listed in SecOJ we assume that Sgr A East accelerates 
nuclei with a power law spectrum with index equal to that of the 7-ray excess 
observed by HESS. 

In SecJU we argue that the spectrum of CRs accelerated in Sgr A East 
may extend beyond the EeV. 

EeV neutrons may be produced as secondary particles of UHE nuclei 
in Sgr A East by two main processes: the Photo-Disintegration (PD) of 
composite nuclei onto the 40 K IR photon background present in that region 
and the pp inelastic scattering of UHE protons onto the dense hydrogen gas. 
In this paper (see Sees El and EJ we estimate the expected neutron flux 
reaching the Earth produced by both processes. We show that PD should 
provide the dominant, if not the unique, contribution. This contribution was 
not computed in a previous work on the subject (7j. 

In SecJ7|we estimate the event rate expected at the Pierre AUGER Ob- 
servatory jTH] due to neutrons from Sgr A East. 

The possibility that EeV neutrons could be produced by a SNR in the 
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GC giving rise to a localised CR excesses around the EeV has been already 
suggested in several papers (see e.g. [T H H5 | ITB|). In particular, Crocker et 
al. [Zj first proposed Sgr A East as a possible high energy 7-ray and EeV 
neutron source, although the contribution of nuclei PD to the neutron flux 
was not determined in their paper. All these works, however, tried to inter- 
pret the EeV CR anisotropy claimed by the AGASA [Hj and SUGAR [T7j 
experiments. This is not what we do in this paper. Actually, in Sec|H] we 
argue that both AGASA and SUGAR anisotropies can hardly be ascribed 
to a neutron emission from the Galactic centre. In particular we show that 
the EeV neutron flux estimated in [Zj is likely to be incompatible with the 
energetic of Sgr A East as constrained by Chandra [5] and HESS [12J ob- 
servations. Here we only use HESS observations to normalise the primary 
spectrum of nuclei accelerated in Sgr A East assuming that this spectrum 
remain unchanged up to UHEs. Although our prediction for the EeV neutron 
flux is much lower than than estimated in [Z] we show that such a flux should 
however be detectable by AUGER. In SecOwe summarise our conclusions. 

2 Physical conditions in the Sgr A East re- 
gion 

Sgr A East has been observed at several wavelengths. Its distance has been 
determined to be d ~ 8 kpc [TSj. 

In the radio this object appears in the form of a non-thermal shell of ellip- 
soidal shape. The major axis of the ellipse has a length of about 10 parsecs 
and runs almost parallel to the galactic plane. The East side of the shell is 
expanding through the dense molecular cloud M-0.02-0.07 . The compressed 
dust /molecular ridge produced by the interaction of these objects give rise 
to a OH maser emission . Yusuf-Zadeh et al. [3 j found some evidence of Zee- 
man splitting in the spectrum of this emission and inferred a magnetic field 
strength of 2-4 mG. Even if such a high field strength is not representative 
of the entire cloud, the synchrotron loss time of relativistic electrons must be 
considerably smaller than Sgr A East age implying that particle acceleration 
is still active in that region. 

Detailed Far InfraRed (FIR) and Medium InfraRed (MIR) observations 
of the GC region have also been performed allowing to identify several diffuse 
and point-like sources. Here we are mainly interested in the diffuse radia- 
tion. In the Nuclear Bulge (R < 300 pc) the IR background appears to be 
dominated by cold dust emission with temperature T ~ 20 — 50 K (FIR). 
This is the true black-body emission as, in the few central parsecs of the 
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Galaxy, and especially in the dense Sgr A East dust ridge, the optical depth 
is much smaller than unity so that the intense UV radiation from the hot 
stars present in that region is effectively thermalised. Warm dust emission 
with temperature up ~ 200 K dominates only in the central parsec of the 
Galaxy. According to the work of Davidson ^H] T ~ 40 K for R ~ 8 pc 
and T ~ 50 K for R ~ 1 pc (see also fUJ E]). Therefore 40 K must be the 
temperature in the NE side of Sgr A East shell. This is also the temperature 
which, according to Gordon at al. 20], is typical for the molecular clouds in 
the Nuclear Bulge. 

The spectacular Chandra [Oj X-ray observations, in the energy range 2 — 
10 keV, with an angular resolution better than 1", confirmed that Sgr A East 
is the remnant of a single type II SN with age between 10 4 and 10 5 years. The 
total luminosity (~ 8 x 10 34 ergs s _1 in the Chandra energy range) and the 
spectrum of this source allowed to establish that its progenitor was a main 
sequence star of mass of 13 — 20 M Q with a metallicity four times larger than 
that of the Sun. The energy of the X-ray emitting plasma was estimated 
to be 10 49 ergs, which is compatible with a total kinetic energy released in 
the explosion of 10 51 ergs, as for a standard type II SN. The X-ray emission 
is concentrated in the central ~ 2 pc of the remnant while no emission has 
been seen coming from the radio shell. 

These observations are compatible with the theoretical predictions for 
a typical SNR of that age which is expected to be in an advanced cooling 
phase. Since the non-thermal radio shell still exists, its present velocity must 
be faster than the sound velocity in the ambient material. According to 
the authors of [S] this implies that the temperature of the ionised gas halo 
surrounding Sgr A East is substantially smaller than 70 K which is consistent 
with the IR observations. The density of the gas halo is about 10 3 cm -3 and 
it is almost homogenous barring the molecular cloud M-0.02-0.07 where it 
becomes as large as 10 5 cm~ 3 [21J. This dense material is swept up by the 
expanding shell giving rise to a thick snowplough in the form of a dust ridge 
surrounding the shell. The dust ridge is substantially denser at the eastern 
edge, where Sgr A East is interacting with the cloud, giving rise to detectable 
molecular line emission (molecular ridge) |21j . 

Soft X-ray and UV and optical radiation are effectively absorbed in such 
a dense medium. Even outside this opaque ridge, the flux of UV and X-ray 
photon density in the surrounding of the radio shell is much smaller than the 
IR radiation density. Indeed, the density of UV photons required to explain 
the gas ionisation in the GC circumnuclear disk has been estimated to be 
< 10 4 cm -3 [22] • An even smaller photon density has to be expected over 
most of Sgr A East region. The X-ray luminosity of Sgr A East observed 
by Chandra implies a negligible photon density over most of such SNR. The 
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ionising source of the gas halo surrounding the GC was not identified by 
Chandra. According to the authors of jS] this source may be related with a 
past activity of the super-massive black-hole in Sgr A* being presently in a 
quiescent phase. The milliGauss field strength measured by Yusuf-Zadeh et 
al. |3j in the Eastern side of Sgr A East may be explained by the adiabatic 
compression of the pre-existing field with an intensity B GC ~ 10 — 100 fiG 
. Polarimetric observations [2~3*j . which mapped the projected magnetic field 
direction in the Sgr A region, support this scenario. Therefore we think it 
is reasonable to assume a magnetic field strength of few milliGauss in the 
molecular ridge crossing M-0.02-0.07, and an order of magnitude smaller field 
in the rest of Sgr A East non-thermal shell. The field coherence length should 
be several parsecs in both cases. 

For what concerns particle acceleration, the weaker field present in the 
latter region is compensated by the larger velocity of the western shell. 
Indeed, the theoretically predicted shock velocity in the gas halo is v s ~ 
6 x 10 7 (t/5 x 10 4 y r ) -0 ' 6 cms" 1 2 J while it is expected to be few times 
smaller in the 10 5 cm -3 dense molecular cloud. Barring the effect of energy 
losses, the maximal CR acceleration energy should be, therefore, roughly the 
same over most of the Sgr A East radio shell. However, hadronic scattering 
in the high density cloud should prevent nuclei acceleration up to the EeV 
by taking place in that region. Acceleration up to the EeV may be also pre- 
vented in the South- West arc, less than 1 pc away from the GC, owing to 
pair production and nuclei PD energy losses (see below) onto the ~ 100 K 
temperature IR background present in that region. 

3 Gamma rays from the GC 

Several experiments observed a 7-ray excess in the direction of the GC. 

EGRET satellite j21| detected an excess centred on the GC with an uncer- 
tainty of 0.2°, corresponding to a projected lenght of 30 pc at the GC distance. 
In the last few years several Cherenkov telescopes allowed to explore the GC 
emission at higher energies. WHIPPLE and CANGAROO II [H] tele- 
scopes, also found a 7-ray excess in direction of the GC in the energy range 
0.2 — 1 TeV. However, the angular resolution reached by these instruments 
was not much better than EGRET. A breakthrough came only recently with 
the HESS telescope which, thank to its stereo-graphic shower reconstruction 
technique, succeeded observing the GC emission up to 10 TeV with an an- 
gular resolution of 1' [T2]. The background-subtracted differential flux mea- 
sured by HESS is F(E) ~ (2.76 ±0.33) x 10" 15 (£/1 TeV) -2 ' 2 cnr^GeV -1 
in the energy range 165 GeV — 6 TeV jT2j. The angular distribution of the 
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signal is compatible with a source of size < 7 pc (95% C.L.), at the distance 
of the GC, centred on a point that, within errors, may coincide either with 
Sgr A* or with the centre of Sgr A East remnant. The spectra observed 
by EGRET and HESS can hardly be connected so that the sources respon- 
sible for the two signals are probably different. Since the excess observed 
by EGRET is within the field of view of HESS but is not observed by that 
telescope, the source responsible for the EGRET signal (3EG J1746-2851) 
must be very faint above 100 GeV. For these reasons we do not consider 
3EG J1746-2851 as a probable site for the acceleration of CRs up to ultra 
high energies. 

Several sources have been proposed to be responsible for the excess ob- 
served by HESS. Among them there are the super-massive black hole in Sgr 
A*, in the galactic dynamical centre, a clumped halo of annihilating super- 
symmetric dark matter in the GC and Sgr A East. In our opinion, Sgr A 
East is the most plausible source for the following reasons: a) the spectrum 
of the GC emission observed by HESS is very similar to that of several other 
SNRs as measured by HESS itself and by other Cherenkov telescopes; b) no 
evidence of a high-energy cutoff was found up to 6 TeV which disfavours the 
dark matter related interpretation of the HESS signal; c) the 7-ray angular 
distribution detected by HESS is slightly Eastwards elongated (see Fig.l in 
[12J); d) the 7-ray flux is compatible with the theoretical estimate based on 
the assumption that such radiation is a secondary product of hadronic cosmic 
rays accelerated in a SNR with the same properties of Sgr A East. 

Due to the relevance of the last issue for the aims of this paper we sum- 
marise here the main passages of such computation. 

Similarly to what done in [2Z1 U\ for 3EG J1746-2851, we assume that 
the high energy 7-rays detected by HESS from the direction of the GC are 
generated by the hadronic scattering of high energy nuclei (mainly protons) 
accelerated by Sgr A East onto the ambient gas nuclei (mainly hydrogen). 
The relevant process is pp pp + n 7T ir ±, ° (n n is the pion multiplicity) with 
the neutral pions decaying into photons. At high energies the differential 
pp cross section can be approximated by a scaling form [2H] and the photon 
emissivity can be written [2HJ 123 Hj 

Q 7 (£ 7 ) ~ n p (E 7 )a™(E )c n H Y 7 (a) , (1) 

where Y 1 (a) is the so-called 'spectrum- weighted moment', or 'yield', of sec- 
ondary photons. In our case Y 7 (2.2) = 0.103 [23 Hj- Reasonably, the mean 
gas density in the Sgr A East region will not differ too much from that of 
the ionised gas halo surrounding the GC which, according to Chandra ob- 
servations, is ~ 10 3 cm -3 (while the gas density is larger than this value in 
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the molecular/dust ridge, behind the shell it is smaller). We assume equipar- 
tition between the total energy of the relativistic particles accelerated by 
Sgr A East and the thermal energy of the X-ray emitting plasma which is 
amounting to 10 49 ergs [9 . Then, by assuming a power law proton spectrum, 
f E\~ a 

Np(E) = Nq I — , and imposing the normalisation condition 
\ E oJ 

f°° 1 

E RC = / dE EN p (E) = -N E 2 . (2) 

Je a - 2 

we find that, for a = 2.2 and E = lGeV, the expected 7-ray flux reaching 
the Earth is 

-15 ( E l \ 2 2 ( n H \ ( E C R \ ! _ 2 _i 



FJE^ ~ 5 x HT i& : n - Z> GeV-Vm-'s 

7V 7; \lTeVj VlO 3 cm- 3 / V!0 49 ergs/ 

(3) 

This is less than a factor two larger than the flux measured by HESS 
which, by accounting for the large uncertainties in the quantities involved 
in this computation, amounts to an excellent agreement between the model 
predictions and the observations. 

We conclude this section by observing that TeV 7-ray attenuation due 
to pair-production scattering onto the IR photon background in the Sgr A 
East region is negligible. The mean IR photon energy is, in fact, 3kT ~ 

1CT 2 ( — — j eV, implying a 7-ray threshold energy of about 20 TeV. The 

effect may be significant only in that portion of Sgr A East shell which is 
within 1 pc from the GC where T ~ 100 K. Gamma-ray scattering onto the 
UV and X-ray backgrounds present in that region was also showed to give 
rise to a negligible attenuation [Zj. 



4 Maximal energies of nuclei 
4.1 Acceleration time scale 

In the first order Fermi acceleration process, the acceleration time scale co- 
incides with the mean time taken by a particle to cycle through the shock 
[2]. This is roughly given by 

t^E) * *m ( 4) 

c u s 

where D(E) is the diffusion coefficient and us is the shock velocity. When 
the particle energy is very high it is reasonable to assume Bohm like diffusion 
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giving the minimum value of the diffusion coefficient D(E) = ^cRl(E). The 
Larmor radius of a nucleus with charge Z is 

R ^ = zTb ~ 1 10 (lo^v) (io^g) ' pc • (5) 

An important condition to be fulfilled is that RL(E max ) does not exceed the 
size of the SNR. 

Using the previous expressions we find 

Nuclei acceleration up to UHEs might be prevented by one of the fol- 
lowing processes: a) hadronic collisions onto the gas particles; b) electron- 
positron pair-production by the interaction with low energy photons; c) 
photo-disintegration onto background photons. The relative relevance of 
these processes depends on the physical conditions in the acceleration re- 
gion. 



4.2 Energy losses due to hadronic collisions 

The gas density in the GC region is considerable larger than the mean ISM 
density. Chandra observations [§] allowed to infer that Sgr A East shell is 
expanding through a ionized gas halo with mean density nn ~ 10 3 cm -3 . 
The hydrogen density can be as high as ~ 10 5 cm -3 in the molecular cloud 
M-0.02-0.07 . The nucleus-proton scattering time scale is 



^hadr 



I A- 2 / 3 10 12 ( UH \ 1 s (7) 

A 2 / 3 cn H Ka pp ~ VlO 3 cnrV ' { ) 



where the pp cross section <j pp (E) and the inelasticity K are given in the 
Appendix. The A 2 / 3 factor in (JJJ) is due to the grown of the effective surface 
taking part to the pA interaction respect to the pp scattering. We can safely 
neglect here the weak energy dependence of a pp and K. A more detailed 
treatment will be performed while determining the spectrum of the secondary 
particles produced by the pp scattering (see SecEJ). 



4.3 Pair-production energy losses 

As we discussed in SecElthe region where Sgr A East is located is pervaded, 
besides by the hydrogen gas, by a background of black-body radiation at a 
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temperature ~ 40 K. Therefore, the photon density in that region is consid- 
erably larger than in the interstellar and intergalactic media where the CMB 
is the dominant radiation background. It is well known in cosmic ray physics 
that a photon background give rise to significant energy losses of UHECRs. 
Besides the PD, which in our case is a useful process and will be consid- 
ered in Sec0 the most relevant loss process for nuclei is the pair production: 
((A, Z)j — > (A, Z) + e + e~) . This process has been studied in details in two 
seminal papers by Blumenthal and Chodorowsky et al. jHI]. Applying 
the results found in those papers we find 



i 2ari Z l m e , , „ . m , 
=^T^ 3 /(7T) 

* l ( r V . , 



-rfilT) — - 1.4 x 10~ i5 s" 
A J w \ 40 K ' 



where 7 = E/Attin- The function 



7/cT 



/(tt)= t^t de0(O-=d — > (9) 



where 0(£) is defined in [3T], is plotted in Fig. [T] 




yT (MeV) 



Figure 1: The function f(jkT) is represented. 

Therefore, it follows from Eq.s (jHJ) and (0) that in Sgr A East pair- 
production losses of ultra-relativistic nuclei are subdominant respect to hadronic 
losses. Photo-pion production also does not play any significant role in the 
energy range of interest for this work. 
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4.4 Maximal energy 



An order of magnitude estimate of the maximal energy reachable by diffusive 
shock acceleration can be determined by imposing the condition 



tacc(E) 



min{t loss (£;), t SNR } 



(10) 



where t acc (E) : t\ oss (E) and tsNR are respectively the characteristic accelera- 
tion time, energy loss time and the age of the SNR. As we showed in the 
above, energy losses in the Sgr A East region are dominated by hadronic 
scattering with a characteristic time given by Eq. ( |7| ). Sgr A East age is 



tsgrAEast ~ 5 x 10 4 yr ~ 10 12 s . (11) 



which, especially for heavy nuclei, is larger than thadr- 

Therefore by substituting © and (J7J) into (|10|) we find that the maximal 
energy reachable by shock acceleration in Sgr A East is roughly given by 



E ~ JL ( ?sl£\ ( B \ 10 i8 eV m 

max " A 2 /3 ^ io-3 ) \ 100 ) [ ] 



We conclude that the acceleration of nuclei up to the ankle of the CR spec- 
trum in Sgr A East is plausible. According to Croker et al, even larger 
energies may be reached in that object if the magnetic field has a component 
non-perpendicular to the shock. 

5 Neutron production by nuclei photo-disintegration 
onto the IR background 

To determine the production rate of neutrons by the photo-disintegration 
(PD) of high energy nuclei interacting with a heat bath of photons we apply 
the results of Puget, Stecker and Bredekamp [32] • Following what done in 
|32j . the photo-disintegration rate for a nucleus of atomic number A (we 
consider only stable nuclear species) releasing i nucleons can be written as 




(13) 



where 7 is the Lorentz factor of the nucleus, and e, e' the photon energy in 
the observer and nucleus rest frames respectively. In general, the PD cross 
section is well approximated by the sum of two terms: a term describing one- 
and two-nucleon PD, which is dominated by the giant dipole resonance, and 
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a non-resonant term which is relevant for the multi-nucleon emission only. 
For 7/cT <C 30 MeV, which is the case in the present context, the latter term 
gives a negligible contribution. The function 

2 

a A ,i{e') ~ ^0(e' - 6^)0(30 MeV - e') W^^E* exp 

i=l 

(14) 

was found to provide a good fit of the experimental data [32j- The definitions 
of the parameters e* i; e^, A^j and £4^ are given in [321 EH]- In TabUJwe 
report their values for the PD of most abundant nuclides with the emission 
of only one nucleon. For the threshold energies we use the values given 
in Stecker and Salamon [SB] which we also report in the same table. The 
strength of the PD reaction, i.e., the integral of cross section over the relevant 

(A - Z)Z _ 2fi 2 

energy interval, is E^ = 5.98 x 10 MeV cm . The normalisation 

J\ 

constant is given by 



Nuclear species 


A 


Z 


e l 


4i 


Ai 


6 


Fe 


56 


26 


11.2 


18 


8 


0.98 





16 


8 


15.7 


24 


9 


0.83 


C 


12 


6 


18.7 


23 


6 


0.76 


He 


4 


2 


20.6 


27 


12 


0.47 




erf 



30 MeV - e' 



0, A,i 



A^/v/2 



+ erf 



z A, 



— e 



th 
A.i 



&A,i/V2 



(15) 



Table 1: PD cross section parameters for the most abundant nuclides as given 
in E21IHH1. The e's and A are in MeV. 



We find convenient to express the photo-disintegration rate in the form: 



rj-\ \ 3 



RaAiT) ^RoUil^p) ®A,i(jT) (16) 



where 



Rq = — (5.98 x 10" 26 cm 2 ) (k 40 K) 3 ~ 5 x 10~ 10 s" 1 . (17) 
2n 2 n- i c 2 

The function & A ,i is defined by 

Qjvix) = (WA) A \ * fJ- x- 3 / de (e-> - l)" 1 J A>i {e) , (18) 

A J*% 

12 



where 

J Ai (e) = dy yexp |-2 ~ J . (19) 

We verified that further contributions to the PD rate which were considered 
in |32j are subdominat in our case. 




1 10 1 10 

ykT (MeV) ykT (MeV) 



Figure 2: The function $A,i{lkT) is represented for 4 He, 12 C, 16 and 56 Fe. 

In Fig. El we represent the function &a,i for several nuclear species. In- 
terestingly, we found that for all relevant nuclei the photo-disintegration rate 
is peaked at jkT « 4 MeV. We will see that this property has relevant 
consequences for the spectrum of the secondary neutrons. In terms of the 
nucleus energy the peak condition reads 

E^ k ~ A x 10 18 1 eV . (20) 

The maximal value of $>A,i is different for different nuclear species. We con- 
sider here only those species which are expected to have the largest abun- 
dance in Sgr A East. The reader can see from Fig|2]that, among the most 
abundant elements, iron ( 56 Fe) nuclei have the highest maximal value of the 
PD rate. However, due to grow of energy losses with A (see Eq|7J), we do not 
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expect heavy nuclei to reach E^ 3, . Furthermore, even disregarding hadronic 
energy losses, a too rapid PD (note that the maximal PD rate of 56 Fe is much 
larger than the inverse age of Sgr A East) would itself prevent Iron nuclei 
reaching E v A eak . This is not the case for 4 He nuclei having a PD peak rate 
which is almost coincident with the inverse of the Sgr A East age. We assume 
here that the physical conditions in Sgr A East (namely the local magnetic 
field intensity and geometry nearby the shock) allow the acceleration of 4 He 
nuclei up to few EeV, corresponding to a value of $4 of (2 3) x 1CT 2 . As 
we showed in Sec. EJsuch conditions are not unrealistic. 

In order to evaluate the flux of secondary neutrons reaching the Earth 
we need to estimate the relative abundance of the main nuclear species in 
the SNR shell. Since we are assuming that all nuclei are accelerated by 
the same mechanism under the same physical conditions, we expect them to 

share the same power-law spectrum, namely Na{E) = f A No > where 

Ja parametrizes the relative abundance of nuclear species in the SNR shell. 
Although, as we discussed in SecJH the spectra of different nuclear species 
may have different UV cut-off energies, for a > 2 this has no significant 
effect on the CR energetic. N Q can be readily determined by imposing the 
normalisation condition (J2J). 

Assuming a uniform cosmic ray density in the acceleration region, the 
neutron emissivity is therefore 

1 2 ( AE \ ~ a 

Q(En) = ^ EX>'M« - 2 )^R [~^J min {R A , h t^ R } (21) 

where = (m p + m n )/2 and V is the volume of the emitting region. We 
assumed a probability 1/2 for the expelled nucleon to be a neutron. We 
used here that any neutron expelled during the photo-disintegration carries, 
in average, a fraction 1/A of the parent nucleus energy. The two nucleon 
emission process is subdominant and we neglegt it in the following. The 
minimum condition in ()21)1 is required since for Ra(^T) > tg^ R any further 
acceleration is effectively inhibited by the PD energy losses. 

According to Chandra observations, Sgr A East metallicity is about four 
time that of the Sun. As a consequence we estimate the iron/hydrogen 
relative abundance to be fee ~ 1CT 3 . The PD condition (i?(56) ~ ^snr) is 
reached for 7/cT ~ 0.7 MeV corresponding to a neutron energy ~ 2 x 10 17 eV. 
Carbon and Oxigen should have a comparable abundance and a smaller peak 
rate so that their PD gives rise to neutrons of energy ~ few x 10 17 eV. These 
neutrons however are hardly detectable due to their smaller time-life and the 
small abundance of the corresponding primary nuclei. The most significant 
contribution to the neutron flux reaching the Earth should come from the 
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PD of 4 He nuclei. Indeed f± ~ 0.1 in Sun, and it is presumably even larger in 
the remnant of a type II SN (unfortunately 4 He is quite difficult to observe). 
Furthermore, as we showed in the above, for the 4 He the PD condition is 
reached just at the energy where the PD rate takes its maximal value so 
that, as follows from Eq. fl2()L the flux of secondary neutrons will peak at 
the energy 

£peak _A_ = x 10 18 / y ^ 

This is just the energy required for the neutron flux to be not suppressed 
significantly by the decay on the way from the GC to the Earth. 

Using (j2"Tj) and accounting for the finite lifetime of neutrons we finally 
find that for a = 2.2 the expected neutron flux reaching the Earth from Sgr 
A East is 

dF(E n ) 27 / U\ ( RoZA*A(E n ,T) \ ( Eqk \ ( d x 2 

dE n X VO-l/V 10~ 12 s- 1 J UO 49 ergs/ V8 kpc 

(23) 

Note that, since the neutron mean time-life is r n = 887 s jSl], the neutron sur- 

(d TTl 
Vl. 

0.43. 

The spectrum of neutrons produced by the PD of 4 He nuclei is repre- 
sented in FiglBl While below 1 EeV the flux is exponential suppressed due 
to the combined effect of the PD threshold and the neutron decay, above 
that energy the neutron spectrum is a power-law of index ~ 5 as an effect 
of the convolution of the primary nuclei spectrum with the E~ 3 high energy 
behaviour of the PD rate. The value f± = 0.1 was assumed drawing FiglBl 
This is a quite conservative assumption, as the relative 4 He abundance in Sgr 
A East is likely to be few times larger. 

We conclude this section by observing that PD of 4 He nuclei onto the 
UV and X-ray backgrounds in the surrounding of Sgr A East cannot prevent 
their acceleration up to UHEs. In fact, in SecElwe argued that the number 
density of photons in that energy range is nuv — n x ^ 10 4 cm -3 . This is 
smaller than the photon density of the IR background at 40 K by more than 
two orders of magnitude. We showed in the above, that the PD rate does not 
depend on the absolute value of the nucleus energy but only on the product 
of its Lorentz factor with the photon energy. As a consequence, the maximal 
PD rate due to UV and X-ray is the same of (fTKj) but for the re-scaling factor 
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Figure 3: The differential flux reaching the Earth of neutrons produced by 
the photo-disintegration of 4 He nuclei in Sgr A East 

nuv, x/^ir — 1CT 2 . In the case of 4 He nuclei this quantity is smaller than 
Sgr A East age by the same amount. 

6 Neutrons from pp collisions 

Besides the PD process discussed in the previous section, EeV neutrons might 
be also be produced by the collision of UHE nuclei with the dense gas in the 
GC environment. The relevant neutron production channel is the charge ex- 
change pp inelastic collision pp — > pn + n n ir, where n w is the pion multiplicity. 
The neutron emissivity is 



tity is experimentally undetermined for E p larger than few hundred GeVs. 
We handle this problem by combining Monte Carlo simulations with an an- 
alytical calculation based on the scaling properties of the differential cross 
section at high energy. In the Appendix we find 




(24) 



where 



d(Tpp\E n , E p ) 
dE n 



is the pp differential inelastic cross section. This quan- 




(25) 
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where a = a pp (E ) and the functions /1(a) and 12(a) are given in (JHHj) . 

Respect to a similar approach followed by the authors in our scaling 
function provides a better fit of the MC data at high values of E n /E p and 
takes properly into account the grown of the cross section at very high en- 
ergies. As we did in the previous sections we assume here that the energy 
spectrum of protons accelerated by the SNR is a power law of index a. 

Then, by assuming a = 2.2 and accounting for the finite lifetime of the 
neutron, we find that the differential neutron flux due to pp scattering which 
should reach the Earth is: 

K(En) H 8xl0 ->(_|^)-"(_j*_) (_*|_) 
(^)M-!^) GeV ~ wvv (26) 

This is represented in Fig0J 




Energy (GeV) 

Figure 4: The differential flux reaching the Earth of neutrons produced by 
pp collisions. 

A comment is in order here. In the appendix we showed that the mean 
value of the neutron elasticity ((x) = (E n /E p )) is ~ 0.05. Therefore, the 
maximal proton energy must be larger than 2 x 10 19 eV to give rise to un- 
suppressed neutron flux. In the derivation in the above we didn't account 
for a UV cut-off in the proton energy spectrum implicitly assuming that this 
is far above the EeV. This assumption, however, is probably unrealistic even 
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for Sgr A East (see SecHJ). As a consequence, a suppression of the neutron 
flux due to pp scattering has to be expected respect to the estimate given in 
126(1 unless Sgr A East is an extremely energetic proton accelerator. 



7 Expected signal at AUGER 



The Pierre AUGER Observatory is an Extensive Air Shower detector under 
construction in Argentina [13]. When completed, it will consist of 1600 Water 
Cherenkov elementary Detectors (WCDs) covering an area of about 3000 km 2 
overlooked by a set of four florescence telescopes. An Engineering Array of 
WCDs has been in operation since 2002 and the detection of several showers 
with energy of the primary ~ 10 18 eV has already announced at the ICRC 
2003 [35]- The AUGER acceptance will be AQ « 4700 km 2 sr, for a zenith 
angle 9 < 45°. The array trigger efficiency at the EeV has been estimated to 
be r] ~ 0.25 !33]. 

The geographic position of AUGER is ideal for observing a UHECR excess 
in the direction of the GC. Since the angular size of Sgr A East is few primes, 
the neutron emission from this object should give rise to a point-like signal 
in AUGER. 1 More precisely, the expected signal should be seen under the 
form of a CR excess in the angular bin corresponding to the position of the 
source. The angular bin covers a solid angle u = 7r5 2 , where 5 is the angular 
resolution of the detector. The AUGER angular resolution 5 is expected to 
be around 1° at 10 19 — 10 20 eV (implying to ~ 1 x 10 -3 sr) but it 5 will 
probably be 2° — 3° (in the final configuration) around the EeV ^3] • 

The isotropic background of ordinary CRs gives rise to an unavoidable 
noise. The CR spectrum in the energy range 4 x 10 17 < E < 6.3 x 10 18 eV 

is [3H] 



Therefore the maximal (at 95%C.L.) differential flux around 10 18 eV incident 
onto the expected AUGER angular bin is 

F-(E) * 3 x 10- ( I5 |^) " (A) 2 GeV-WV 1 . (28) 

It is comforting that the expected neutron flux from Sgr A East produced 
by the photo-disintegration of 4 He nuclei (see Eq. |2*3*|) is comparable to the 

1 In principle, EeV protons produced nearby the Earth by neutron decay may give rise 
to a smeared signal superimposed to the point-like excess |36l I37| . This secondary signal, 
however, is probably undetectable. 




-3.20±0.05 



J CR (E) = (9.23 ± 0.65) x 10 



(GeV sr cm 2 s) 



-i 



(27) 
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isotropic CR flux around the EeV incident over the AUGER angular bin. 
As a consequence, AUGER may observe a significant UHECR excess in the 
angular bin corresponding to the position of the GC. 

As a benchmark, we give here the expected rates of showers incident over 
an angular bin of 2.5°, in the energy range 1 x 10 18 <E<3x 10 18 eV, due 
to the PD and pp neutrons respectively (accounting for the AUGER trigger 
efficiency) 

A> PD ~ 132 yr- 1 N pp ~ 72 yr" 1 , (29) 

A^pd has been computed assuming the 4 He relative abundance f± = 0.2 (see 
the discussion at the end of SecEJ). These rates have to be compared with a 
CR background rate in the same angular and energy bins which is 

A> CR ~ 440 yr" 1 . (30) 

Since the expected signal is a sizeable fraction of the background a significant 
excess should be detectable in the direction of the GC. 

We conclude that, under the conditions that we assumed in this paper, 
the signal due to EeV neutrons produced by Sgr A East should not be missed 
by AUGER. 

8 Discussion 

The possibility, discussed in this paper, of testing if nuclei acceleration takes 
place in Sgr A East up to Ultra High Energies (E > 10 18 eV) by looking for 
EeV neutrons coming from the GC direction, was also considered in a recent 
paper by Crocker et al. [7j. Our work, however, differs from that of Crocker 
et al. in a number of relevant aspects. 

First of all, the aim of Crocker et al. was to interpret the anisotropy in 
the CR angular distribution observed, around the EeV, by the AGASA [Hj 
and the SUGAR jTTj experiments in a direction relatively close to that of the 
GC. This is not our attitude. Indeed, even accepting AGASA and SUGAR 
results as evidences of real CR anisotropies, which is tricky due to their rather 
low statistical significance, it is, in our opinion, rather hazarding to identify 
Sgr A East as the source responsible for these signals. Indeed, neither the 
AGASA nor the SUGAR UHECR excesses are in the direction of the very 
GC. The GC is not even in the "field of view" of AGASA being below its 
geographic horizon. 2 While SUGAR is in the right geographic position for 

2 We have to note here that an extended signal due to the protons produced by neutron 
decay has to be expected [23 02] ■ Although, this signal might explain the anisotropy 
observed by AGASA in terms of a neutron source located in the GC, the required neutron 
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observing the GC, no excess was observed by this experiment from the true 
centre of the Galaxy. The SUGAR peak was 7.5° away from the GC and it 
does not even appear to be centred on the galactic plane. This offset is too 
large to be due to a possible systematic pointing error of the experiment. For 
these reasons we think that Sgr A East cannot be responsible for the AGASA 
and SUGAR anisotropies. It is important to notice that SUGAR was not 
sensitive enough to be able to detect the neutron flux which we predicted 
in Sees El and El to come from Sgr A East. Therefore, the lack of signal 
in the direction of the GC in the SUGAR data set does not contradict the 
hypothesis raised in this paper. 

The normalisation of the relativistic particle spectrum in Sgr A East 
which was adopted in [7j was based on a combined fit of the low energy 
7-ray signal observed by EGRET [21] in the direction of the GC and the 
AGAS A/SUGAR anisotropy data. As we explained in SecJHl and noted 
also by the authors of [Zj, the 7-ray spectrum of the source responsible for 
the EGRET excess cannot extend up to very high energies. It is, there- 
fore, quite unnatural to assume that this source accelerates CRs up to the 
EeV. Although in (7j the authors also (marginally) consider a HESS + 
AG AS A/SUGAR fit, in order this to work they need to adopt a power index 
1.97 for the proton spectrum which is almost 3a lower than the HESS best 
fit [12] (as we showed in SecJHl the spectrum of high energy secondary pho- 
tons is expected to have the same slope as that of primary protons). Such a 
discrepancy may become more serious if new data from HESS would prolong 
the GC 7-ray spectrum observed so far to even higher energies. 

Last but not least, the contribution of nuclei photo-disintegration to Sgr 
A East neutron emission was not determined in [Zj. We showed, however, 
that this process should provide the dominant, if not the unique, channel for 
the production of EeV neutrons. 

In this paper we did not discuss the contribution to the synchrotron radio 
emission and to the inverse Compton (IC) low-energy 7-ray flux from Sgr A 
East which may be produced by secondary leptons. The possible relevance of 
these emissions was studied in j2Zj. The authors of [211 assumed that the 7- 
ray flux observed by EGRET [24J is a by-product of pp inelastic scattering in 
Sgr A East and used EGRET measurements to normalise the primary proton 
flux. On the basis of these assumptions, they estimated the expected radio 
and IC secondary fluxes showing that these emissions are compatible with 
the observations. Since the normalisation of the relativistic nuclei spectrum 
that we used in this paper (which is based on HESS observations) is much 

luminosity of such source, which was determined in [37| . is too high to be compatible with 
its possible identification with Sgr A East. 
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lower, we are confident that secondary emissions are even harmless in our 
case. 

We also did not discuss here the possible effect of post-acceleration diffu- 
sion on the primary nuclei spectrum. In principle, since diffusion is an energy 
dependent process, deviation from a single power-law behaviour may arise 
in the nuclei spectrum. Some possible consequences of high energy proton 
diffusion in the GC region has been recently discussed in Here we im- 
plicitly assumed that above several TeV (the energy required to produce the 
TeV gamma-rays observed by HESS) protons spectrum around Sgr A East 
coincides with the acceleration spectrum. 

9 Conclusions 

In this paper we estimated the flux of neutrons reaching the Earth from Sgr 
A East under the educated hypothesis that nuclei are accelerated beyond 
the EeV in that SNR. We showed that neutrons with the required energy 
(~ lEeV) can be naturally produced by the photo-disintegration (PD) of 
composite nuclei onto the IR background radiation present in the GC re- 
gion and by the pp scattering of UHE protons onto the dense hydrogen gas 
surrounding Sgr A East. We found the former process to be the dominant. 
Amazingly, the PD rate of nuclei onto the 40 K photon background peaks 
just at the energy required to produce EeV neutrons which can reach the 
Earth from the GC before decaying. Due to their largest abundance respect 
to other composite nuclides, 4 He nuclei should give the main contribution 
to the neutron flux. While the energy threshold for pp scattering is much 
lower, due to the low neutron elasticity (the fraction of the primary proton 
energy going to the neutron) the energy required to primary protons in order 
to produce EeV neutrons is, in average, higher than 10 EeV. Such a high 
energy may be hardly reachable even in Sgr A East so that the contribution 
of pp inelastic scattering to the EeV neutron flux may be absent. 

The signal produced by the EeV neutrons from Sgr A East should be 
easily disentangled from the CR isotropic background thanks to its particular 
spectrum and its point-like nature. As we showed in SecEl the neutron 
spectrum should start to rise rapidly above few x 10 17 eV, to peak just below 
the EeV, and become a power law steeper than the CR spectrum at higher 
energy. The reader should take in mind, however, that a UV cut-off may 
come-in somewhere above the EeV as the maximal acceleration energy in Sgr 
A East cannot be indefinitely high. Presumably, only one energy bin around 
the EeV will be interested by the excess. The clearest signature of the neutron 
flux from Sgr A East should be given, however, by the point-like nature of 
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the excess and by its coincidence with the actual position of this SNR in the 
sky coinciding with that of the GC. By assuming that nuclei accelerated in 
Sgr A East have a power law spectrum with exponent —2.2, as suggested 
by HESS gamma-ray observations, and that this spectrum extend steadily 
beyond the EeV, we determined the flux of secondary neutrons expected to 
reach the Earth and showed that the Pierre AUGER Observatory should be 
able to detect it. 

We conclude by observing that several of the results that we derived in 
this work for Sgr A East should also be applicable to some other SNRs. 
As we discussed in the Introduction, the most promising neutron/gamma- 
ray emitters as those SNRs interacting with dense molecular clouds. The 
fast improvements in gamma-ray astronomy and in the extensive air shower 
reconstruction techniques may allow soon a systematic search of galactic 
EeV-trons. 
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Appendix 

The pp total cross-section determined by combining collider and UHECR 
(PP air scattering) data is 



where Z pp = 35.48 x KT 27 cm 2 , B = 0.310 x 10~ 27 cm 2 and si = (5.38 GeV) 2 . 

Clearly, we are interested only in the inelastic pp scattering. The inelas- 
ticity K is defined by a l p n p el = Ka pp From the experimental data of the pp and 
pp scattering K ~ 4/5. For s > 10 GeV 2 there is no evidence of a variation 
of K with energy. 

The neutron emissivity can be written in the form 




(31) 




(32) 
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dP n 

where n n is the neutron multiplicity and ——(E n , E p ) is the normalised prob- 

dE n 

ability distribution to produce a neutron with energy E n . 

Since neither experimental data nor undisputed theoretical calculations 
are available allowing to determine these quantities at UHEs, we resort to 
Monte Carlo simulations. We used the Pythia generator [12] to simulate pp 
inelastic scattering for several values of E p in the range 1 — 10 3 TeV under 
the constraint that the total cross section is that given in (|31j) . We found 
the value of the neutron multiplicity to be n n = 0.24 almost independently 
on E p . The simulation also allowed to determine the neutron energy distri- 
bution. Again, the result was practically energy independent meaning that 

differential cross section — — can be written in a scaling form: 

dE n 

do dP <j mc ^ 

Vinein n —^-(E n ,E p ) = Yl n -g- fl(x) , (33) 



where x = E n /E p . The best fit to the Monte Carlo data gives 

h(x) ~ 0.064 (1 - x) 2 + 0.094x^/T^x . (34) 



This result is similar to that found by Drury et al. pi] but for the last term, 
being absent in their work, which allow a better fit of the MC data for x — > 1. 
The mean value of x, the so called neutron elasticity, is 

(x) = [ dx h(x) ~ 0.05 , (35) 



Before inserting the previous results in (j32j) we rewrite the inelastic cross 
section in the form 



° P f(E p )=o, (l + Dln(|k) ) , (36) 



where we defined a = a pp (E ) ~ K x Z pp , and D = 2B/Z PP = 1.7 x 10 2 
where Eq = 1 GeV. The neutron emissivity can then be written as 



h(a) + I 2 (a)ln 



Ql P (E n ) = cn H n p (E n )a n n K 
where 

h(a)= [ dxx a - 2 h(x){\ -D\xi(x)) 



(37) 



i 

Q-2j 



(38) 



h{a) = D / dxx a ~ 2 h(x) . 
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Clearly, a regularization procedure is needed in order to avoid singularities 
near x = 0. Obviously the correct regularization is made when one considers 
the finite mass of the neutron, which forces the integration domain to be 

m n 
— , 1 
E p 
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